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We re-analyse experiments on a foam sheared in a two-dimensional Couette geometry 
[Debregeas et al, Phys. Rev. Lett. 87, 178305 (2001)]. We characterise the bubble defor- 
mation by a texture tensor. Our measurements are local in time: they show two regimes, 
one transient and one stationary. They provide both the average and fluctuations of the 
anisotropy. Measurements are also local in space: they show that both the deformation 
and the clastic contribution to the stress field do not localise, varying smoothly across the 
shear gap. We can thus describe the foam as a continuous medium with elastic properties. 



1. Introduction 

A liquid foam exhibits "complex" behaviour under stress: it is elastic for small de- 
formation, plastic for large deformation, and flows at large deformation rates (Weaire & 
Hutzler 1999), (Jiang et al. 1999). It is a model for the study of complex fluids: its typical 
length scale permits direct observation. In particular, a foam with only one bubble layer 
(so-called "two-dimensional foam" (Weaire & Hutzler 1999)) is easy to image, and image 
analysis yields information on all the geometrical properties of the foam. 

Recently, Debregeas et al. (Debregeas, Tabuteau & Di Meglio 2001) quasistatically 
sheared a fairly dry foam in a two-dimensional Couette geometry (Fig. 1). They studied 
the velocity field and showed that it localises in a shear band: its changes of value 
happen as a steep gradient on a layer much thinner than the gap, rather than over 
a long distance. It decreases exponentially across the gap (ri — r ), as exp [— (r — r )/d]. 
Here r is the distance to the wheel center, vq the radius of the inner wheel, r\ the 
radius of the outer wheel, and d <C (ri — ro) is a characteristic length comparable to 
one bubble diameter. Such shear-banding has been explained and numerically simulated 
by the combination of frozen disorder, elastic behaviour at low deformation, and local 
discrete bubble rearrangements, leading to an instability more likely to occur at the inner 
side where fluctuations are larger (Kabla & Debregeas 2003). 

What dominates the macroscopic behaviour of the material: averages or fluctuations? 
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Figure 1. Run 1 (transient regime): foam under steady quasistatic shear in a Couette geometry. 
Pictures of the foam correspond to different values of applied shear strain 7: (a) 0, immediately 
after switching from counter-clockwise to clockwise rotation; (b) 0.15; (c) 0.5, when the inner 
boxes have reached the yield strain (see Fig. 5); and (d) 7.6, well into the stationary regime. 
The cross-component r6 of the texture tensor, defined in eq. (3.11), is positive (a), zero (b), 
negative (c), and stationary (d). 



This question is important in the rheology of foams, as well as in other disordered ma- 
terials. Briefly, if the signature of fluctuations averages out at large scales, a "thermo- 
dynamical limit" exists in which the average dominates; a constitutive equation relating 
the stress and strain fields can then describe the material as a continuous medium (like 
a solid or a liquid), ff the fluctuations have the same order of magnitude as the average 
at all scales, they can dominate the physical behaviour, which need not correspond to 
any continuous medium. Whenever the fluctuations dominate, we can ask whether they 
correlate (leading to large stress drops) or have a white spectrum. The answer need not 
be the same for different materials (say, foam and granular materials) nor for different 
quantities (say, velocity and deformation). 

Similarly, the existence and origin of localisation raise different questions: the effect 
of boundary conditions, the respective roles of strain and strain rate, the precise role 
of fluctuations, the differences between various (disordered) materials, the distinction 
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between sharp (Coussot et al. 2002) or exponential (Debregeas, Tabuteau & Di Meglio 
2001) decreases, the dimension of space (2D or 3D), and even the field itself (velocity 
or deformation). We would like to extract additional physical information, beyond the 
velocity field. 

Most suggestions to quantify statistically the anisotropy and the deformation of the 
internal structure of the foam, that is, the bubble wall network, are based on scalar 
quantities (Elias et al. 1999). Since the deformation encompasses information on both 
anisotropy and orientation, it is actually a tensor (Landau & Lifschitz 1986), (Alexander 
1998). Recent suggestions of possible tensorial descriptions based on bubble details are 
promising (Rcinclt & Kraynik 2000), (Kraynik, Reinelt & van Swol 2003), (Ball & Blu- 
mcnfcld 2002), (Blumenfeld 2004). Here, we use the "texture tensor" (eq. 3.11), which 
has the following advantages (Aubouy et al. 2003): (i) it is purely geometrical, indepen- 
dent of stresses and forces; (ii) it applies to both small and large numbers of bubbles; 
(iii) it is local in time and space (Jiang et al. 2000), so it can characterise a sub-region of 
the foam; (iv) it applies to disordered or ordered, 2D or 3D materials, in elastic, plastic 
or fluid regimes. 

The plan of this paper is as follows. Section (2) reviews the original experimental set-up, 
and our image analysis protocol. The theoretical section (3) reviews the definitions and 
measurements of stress, texture and strain tensors; it can be read separately. The next 
section (4) presents our analyses of experiments. The first (section 4.1), resolved in both 
time and space, studies the transient regime before a stationary flow is established. The 
second (section 4.2), resolved in space but averaged in time, shows that in stationary 
flow, both the bubble anisotropy and the elastic contribution to the deviatoric stress 
field vary smoothly across the Couette gap; they do not localise near the inner rotating 
wheel. The third (section 4.3), again resolved in time, shows that in stationary flow the 
fluctuations are larger near the inner wheel. In section (5), we discuss whether we can 
describe the foam as a continuous medium with elastic properties. 



2. Materials and Methods 

2.1. Experiments 

Georges Debregeas generously provided experimental data obtained in two different runs. 
For experimental details see Ref. (Debregeas, Tabuteau & Di Meglio 2001). Briefly, the 
set-up consists of an inner shearing wheel and an outer fixed one, with no-slip conditions 
due to tooth-shaped boundaries (tooth depth: radius 1.2 mm). The foam is confined 
between two transparent plates separated by h = 2 mm spacers (Fig. 2(d)); this foam 
thickness is smaller than the distance between bubble centers (see below), enforcing a 
single layer of bubbles (Cox, Weaire & Vaz 2002). The fluid fraction, measured by the 
weight of water introduced in the cell, is $ = 5.2%. The bubble size dispersity is ~ 30%, 
enough to prevent crystallisation. Coarsening and size sorting under shear are negligible 
during the experiment. A CCD digital camera records the position and the shape of the 
bubbles. 

The inner wheel rotates slowly enough (V w heel — 0.25 mm.s -1 ) that the flow is well 
in quasistatic regime: the foam passes through a succession of equilibrium states. All 
averaged static physical quantities, like the elastic strain or stress, remain the same if 
the applied shear rate V w heei/{ri — r ) varies. Kinetic quantities, like the velocity field, 
the rate of Tls, and the input power and dissipation, scale as the shear rate. This scaling 
is typical of solid-like friction (only at higher velocity does the foam exhibit a fluid-like 
friction, when the input power and dissipation scale like the square of the shear rate). As 
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required for comparison to other experiments, we thus express all results as a function 
of the total applied shear strain: 

Vwheel t /r> -t \ 

7= ■ 2-1 

n - r 

Note that 7 is proportional to time t, not to the local shear in each box. 

Run 1 (unpublished data) contains a transient regime. To prepare the foam, the inner 
disc is rotated counter-clockwise, until a stationary regime is reached. Then, at an arbi- 
trary time, chosen as the origin (t = 7 = 0), the shear direction is switched to clockwise 
(Fig. 1), the experiment begins and 500 pictures are recorded. The r6 components of the 
deformation thus switch from positive to negative (Fig. 3). Each picture shows 200 bub- 
bles, in an angular sector of the experimental cell (Fig. 1). The average distance between 
the centres of mass of neighbouring bubbles is 3.0 mm, with a standard deviation of ±0.5 
mm, corresponding to a mean bubble wall length of (£) = 1.7 mm. The internal radius is 
ro = 71 mm (same for both runs) and the external radius is n = 112 mm (different for 
each run). Thus 7 = 1 corresponds to a wheel displacement of V w heei t = r\—r^ = 41 mm: 
i.e. 15-20 bubble diameters. 

Run 2 (published in Ref. (Debregeas, Tabuteau & Di Meglio 2001)) focusses on the 
stationary regime. The preparatory rotation is clockwise. Then, at an arbitrary time, 
chosen as the origin (t = 7 = 0), it is switched to counter-clockwise (Fig. 2). Pictures 
are recorded only after a full 27r-turn, corresponding to 7 = 8; the rO components of 
the deformation are thus positive (Fig. 6). A quarter of the Couette cell is filmed, with 
700 bubbles in each picture. To improve the statistics, the camera is then displaced 
to record successively the four quarters; 2000 pictures are obtained as four 500-images 
movies. Data are measured and averaged over the four quarters of the foam. The average 
distance between the centres of mass of neighbouring bubbles is 2.4 mm, with a standard 
deviation of ±0.4 mm, corresponding to a mean bubble wall length of {£) = 1.4 mm. 
The internal radius is r = 71 mm (same for both runs) and the external radius is 
r\ = 122 mm (different for each run). Thus 7 = 1 corresponds to a wheel displacement 
of 51 mm: i.e. 20-25 bubbles diameters. 



2.2. Image analysis 

We analyse the images using a home-made extension module to the public software NIH- 
Image. We threshold the image, and skeletise the bubble walls to a one-pixel thickness. 
We identify the points where three walls meet ("vertices"). We replace each wall by the 
vector £ linking its vertices (white lines in Fig. 2a, top right). 

To respect the symmetry of the Couette cell, this wall vector £ (Fig. 2b) is expressed in 
local polar coordinates (r, 9) (Fig. 2a). This neglects the variation of the polar referential 
between both its end points. This assumption is correct as long as £ <C r. 

For the same reason, we choose equi-spaced orthoradial, circular boxes (Fig. 2a, top 
right) to divide the foam and measure the deformation using spatial averages to improve 
the statistics. In the quasistatic regime, the bubble deformation remains very small and 
barely visible by eye. After trying various box sizes, we choose 8 boxes as a compromise 
between the quality of the statistics and the level of visible details. 

For Run 1, each box is 2.7(£) wide, and box number i is at position r\ = ro + [0.4 + 
2.7 (i — 0.5)]{£). The number of bubbles in each box grows linearly, typically from 17 
bubbles in the 1st box, to 32 in the 8th box. 

For Run 2, r\ = r + [2.8 + 3.1(i - 0.5)](£). Due to the shape of the pictures (Fig. 2a), 
the number of bubbles first increases with i (from typically 105 bubbles in the 1st box), 
then decreases (to 41 in the 8th box). 
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Figure 2. Run 2 (stationary flow): image analysis and notations, (a) Bottom left: top view of 
the original experiment by Debregeas et al. (Debregeas, Tabuteau & Di Meglio 2001), rotating 
counter-clockwise. Top right: image analysis. The reconstructed network of vertices linked by 
vectors (white) and the boundaries of the 8 equi-spaced concentric circular boxes used as rep- 
resentative volume elements (RVE) for spatial analysis (black circular arcs) are superimposed 
on the snapshot of the experiment. We exclude the regions near both wheels, where the bubbles 
are not entirely visible, (b) We denote by I the vector linking both vertices of a given wall, (c) 
If one of these vertices is outside of the box (defined in (a), top right), we denote by / £ [0, 1] 
the fraction of the bubble wall length inside the box; / = if the wall is entirely outside, / = 1 
if it is entirely inside, (d) Side-view sketch of a vertical bubble wall cross-section, between two 
horizontal glass plates. 



3. Definitions and measurements of deviatoric elastic stress and 
strain 

3.1. Capillary contribution to the stress 
3.1.1. Theoretical definition 

In fluid foams, the three main contributions to stress a are the capillary term 7j cap 
due to the tension of the liquid walls, the pressure of the gas within each bubble, and 
the viscous dissipation. In what follows, the "deviatoric" part of the capillary stress 
refers to the traceless tensor, i.e. a cap — Tr (a cap ) 1/2, where I is the identity tensor 
in 2 dimensions. It consists in the shear stress a^ p and the normal stress difference 

n-cap _ cap 
°rr °G0 ■ 

In the extreme "wet foam" limit <& ~ 1, bubbles are well separated from each other. 
The capillary contribution in a representative volume element (RVE) V which contains 



6 



E. Janiaud and F. Graner 



N bubbles with surface tension Y is (Roscnkildc 1967), (Batchclor 1970): 



n cap = -j-J- £> (i — n®n\ ds, 



<7 = V ^7 (i&i-wi) ds ) 




(3.1) 



Here N/V is the density of bubbles; S is the surface of each bubble; n is the unit vector 
normal to the surface element ds; and (.)y is the average over bubbles in V . 

In the other extreme, "dry foam" limit $< 1, pairs of bubble surfaces merge into a film 
with surface tension 2T. The capillary contribution is (Khan & Armstrong 1986), (Khan 
& Armstrong 1987), (Rcinclt & Kraynik 2000), (Kraynik, Reinelt & van Swol 2003): 



The correction of order <E> 2 is a contribution of the integral over the surface S of its mean 
curvature (Mecke & Stoyan 2002). Since the mean curvature is mostly concentrated 
in bubble edges ("Plateau borders"), this contribution is physically equivalent to the 
recently identified line tension which affects the bubble edges (Kern & Weaire 2003), 
(Geminard et al. 2004). 

3.1.2. Operational definition in two dimensions 

A true 2D foam has an interesting property: at equilibrium, or in a quasistatic regime, 
the coarse-grained Laplacian of the pressure is zero (Graner et al. 2001), (Graner 2002). 
The small scale fluctuations of bubble pressures, according to bubble side numbers, aver- 
age out to zero at large scales (Weaire & Hutzler 1999). In the present Couette geometry, 
where no large-scale pressure gradient exists, pressure does not vary significantly (Weaire 
et al. 2004). Hence, in-plane curvatures remain much smaller than the inverse of bubble 
sizes. They contribute to the stress through their algebric average (not their second mo- 
ment), which is small along any path within the foam: curvatures are counted positively 
if the path crosses the wall from the concave side (Graner 2002). They thus induce only 
a negligible correction to the deviatoric capillary stress, with respect to straight walls 
(Asipauskas et al. 2003). 

However, adapting the above definition (eq. 3.2) to measure deviatoric capillary stress 
fields in bubble monolayers requires care. In the present case, the curvature radius in the 
third dimension, not visible in images, is smaller or equal to the half distance between 
plates (Fig. 2(d)). Its inverse, namely the curvature in the third dimension, is larger than 
the curvature visible in the image. This case is thus intermediate between two limiting 
bubble monolayers: bubble rafts, and Langmuir foams. 

In bubble rafts, stress has already been measured, in a circular Couette geometry, 
by mechanical measurements of the force at the boundaries (Pratt & Dennin 2003), or 
indirectly by local measurements of bubble ellipticity (Ybert & di Meglio 2002). Since 
bubble top surfaces are curved (they are approximately spherical caps), direct local mea- 
surements of stress are difficult. 

At the other extreme, Langmuir foams are one molecule thick: their aspect ratio is of 
the order of nanometers divided by centimeters, with no well-defined curvature in the 
third dimension. They are thus probably the closest approximation to a true 2D foam (as 
long as 3D viscosity effects remain negligible) , and capillary stress fields can be measured 
as described below (eq. 3.5) (Courty et al. 2003). 

In the present case, for a single layer of bubbles (Cox, Weaire & Vaz 2002) in the 




(3.2) 
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quasistatic regime, the dry part of the film remains vertical (Fig. 2(d)). By symmetry, 
we can assume that the curvatures on both sides of a film are almost opposite, so that 
their effects almost cancel each other. We thus neglect the effect of curvature in the third 
dimension. 

We simplify further eq. (3.2) using two approximations. First, although strictly speak- 
ing it applies to closed surfaces (complete surfaces of bubbles), we extend it to averages 
on open surfaces (walls between bubbles), assuming that any possible non-compensated 
term (due for instance to bubbles at the foam boundary) has a negligible contribution 
to the average (Kraynik, Rcinclt & van Swol 2003). Second, on closed surfaces in 2D, 
the integral of the mean curvature only contributes a constant, so that the O ($ 2 ) term 
need not be considered here; in what follows, on open surfaces, we again neglect the 
contribution of possible non-compensated terms. 

Finally, eq. (3.2) simplifies to an expression which has already been used (but not 
always justified) for various types of bubble monolayers (Weaire & Hutzler 1999), (Kabla 
& Debregeas 2003), (Asipauskas et al. 2003), (Courty et al. 2003): 

* cap = v^ f ( AS ^' (3 ' 3) 

I 

Here (Fig. 2) the 2D RVE V is a surface area; the £s are the chords linking two connected 
vertices, approximating bubble walls; / is the fraction of the wall length contained in V; 
e = £/£ is the unit vector tangent to I. The line tension A is equal to Th multiplied by a 
geometric constant which depends on $ (Fig. 2(d)): it varies from 2 for a dry foam (two 
flat interfaces of length h) to tt for a wet foam (two half-circles of diameter h) . 

3.1.3. Comments on this definition 

Eq. (3.3) is analogous to the expression for the elastic stress in a (2D or 3D) network 
with a density p of links I undergoing two-body interactions with tensions f (Kruyt 
2003): 




In the continuum limit, the average (.)y is the average of the wall probability distri- 
bution function, and this expression for the stress is exactly the same (see for instance 
ref. (Kruyt 2003) or the appendix of ref. (Aubouy et al. 2003)) as the classical definition 
based on the links crossing the boundaries of V (Landau & Lifschitz 1986). This equiv- 
alence validates a posteriori our application of eq. (3.2) to open surfaces. In practice, 
experimental data have a finite size: the number of links used in the calculation, and 
so the signal/noise ratio of the resulting measurement (Asipauskas et al. 2003), is much 
larger if we use all links in V, rather than only the links crossing its boundaries. 

Eq. (3.4) shows that, in a network of points which interact through a pairwise potential, 
we can measure the stress from an image if we know the expression for the two-points 
force law f = f{t) explicitly (for instance Lennard- Jones, or harmonic) . 

For quasistatic 2D foams (eq. 3.3), all tensions have the same modulus A; hence t(£) = 
Ae. The image provides a cap /X, which is enough for our purposes: 
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Here p is the area density of bubble walls: the number of bubble walls in V weighted by 
/, divided by V; (.)v is the average over the walls in V, also weighted by /. 

3.2. Status of the elastic strain tensor 
3.2.1. Macroscopic definition of the elastic strain 

A macroscopic definition of clastic strain (Macosko 1994), (Farahani & Naghdabadi 
2000) is usually a function of the deformation gradient F: 



dx° J 



(vf), 



1 % ( 3 - 6 ) 

where x is the current position of a point initially at position x° . This deformation 
gradient F coincides with the identity / in a solid translation, and has a determinant 1 
in a volume-conserving transformation. 

From F, we can define two symmetric tensors U and V through a rotation R (Phan- 
Thien 2002): 

T = "RU = V% (3.7) 

They are called the right and left stretch tensors. They coincide with the identity I in 
a solid transformation (translation and rotation), and also have a determinant 1 in a 
volume-conserving transformation. 

In turn, they determine a right and a left deformation tensors, called the Finger tensor 
and the Cauchy-Green tensor, respectively (Phan-Thien 2002): 

V = T t f = U 2 . (3.8) 

They too are symmetric, coincide with the identity / in a solid transformation (translation 
and rotation), and have a determinant 1 in a volume-conserving transformation. 

Finally, a strain measure tensor E is for instance equal to (j3 — I^j /2, — I^j /2, or 
more generally to one of the Seth-Hill strain measure tensors based on the successive 
powers of U or V: 

— m zz 

=M U -I 

k = 7 

m 

=m = 
=(m) V - / 

or E = (3.9) 

m 

where m is a positive or negative integer (Farahani & Naghdabadi 2000). They all tend 
towards the linear elasticity definition (the gradient of the displacement field) in the limit 
of infinitesimal transformation. 

The particular case of m = corresponds to the logarithm-based definition (Hoger 
1987) and generalises the Hencky strain tensor (Tanner & Tanner 2003) to arbitrary 
deformations: 

=(0) . == 1 . == 
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This true strain tensor E correctly quantifies both volume-invariant deformations 
(which appear as traceless strain tensors) and dilations (which appear as the logarithm 
of the volume, as, for example, in gases). 

Constitutive equations of a material, that is stress-strain relations, can use either of 
these tensors; some arguments favor right tensors (Farahani & Naghdabadi 2000), (Larson 
1997) while other favor left ones (Macosko 1994). 

3.2.2. Micromechanical definitions of elastic strain 

A microscopic definition of elastic strain (Kruyt & Rothenburg 1996), (Liao et al. 
1997), (Alexander 1998), (Ball & Blumenfeld 2002), (Blumenfeld 2004), (Kruyt 2003), 
(Goldhirsch & Goldenberg 2002), (Zimmerman 1999), (Aubouy et al. 2003) adopts a 
different point of view: it constructs the macroscopic elastic strain from the positions of 
individual "sites" {i.e. atoms, grains, or vertices of a network), just as stress is constructed 
from individual forces, see eq. (3.4). 

For instance, Zimmerman considers elastic deformations of 2D or 3D ordered atomic 
crystals (Zimmerman 1999). He constructs separately the stress (from atomistic poten- 
tials) and the elastic strain (from displacements). He shows that, in the ordered cases he 
considers, his strain yields a continuum-mechanics limit identical to that derived from 
the Cauchy-Green tensor. 

Kruyt (2003) expands the analysis he used in a series of papers with Rothenburg, 
which, for instance, predicted the elastic moduli of disordered or ordered 2D granular 
materials (Kruyt & Rothenburg 2002). He constructs separately the stress from indi- 
vidual forces (obeying equilibrium conditions) and the elastic strain from individual 
displacements (obeying geometric compatibility conditions) small enough to avoid re- 
arrangements (Kruyt 2003). 

Goldhirsch and Goldenberg (2002) rewrite elasticity theory on a microscopic basis, 
take the continuous limit, and show its compatibility with classical elasticity, including 
the energetics and conjugation relations (Goldhirsch & Goldenberg 2002). In particular, 
they show that their definition of elastic strain is better than those based, for instance, 
on fits of the displacements to a macroscopic strain field (Liao et al. 1997) and that, even 
for granular materials, microscopic fluctuations progressively average out when going to 
macroscopic scales. 

3.3. Extension to plastic and fluid regimes 

Macroscopic and micromechanic elastic strain tensors become difficult to define after the 
first topological rearrangement. They finally lose their meaning after a few rotations of 
the Couette-cell inner wheel, and mathematical singularities appear. This breakdown of 
formalism is counterintuitive; experimentally, the mechanical behaviors of elasto-plasto- 
viscous materials do not appear to change discountinuously. 

We need to treat the (even large) elastic deformation, independently from the flow 
involving rearrangements of the connections between neighbouring sites (changes in net- 
work topology). The total applied shear strain decomposes into plastic and elastic strains. 
The plastic strain is much larger, since it increases linearly with time (its derivative, cor- 
responding to the velocity gradient, is constant in steady flow). The elastic strain is much 
smaller: here, in steady flow, it remains constant; but it is physically important, because 
it characterises the current state of the foam. 

We use below the approach of Aubouy et al. (2003), which consists in three steps. It 

first characterises the current state of the material by a "texture" (or "fabric") tensor M, 
built from statistical averages over microscopic positions, which exists in all regimes. It 
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then defines a reference state characterised by a reference tensor Mq. It finally constructs 
the elastic strain by a comparison of both tensors. We now review these definitions. 

3.3.1. Texture tensor 

The pattern of the walls in a region of the foam can be statistically characterised by 
the local texture tensor M (Aubouy et al. 2003). In 2D, it writes (again with averages 
weighted by /): 

M ij = (i i t j )v = ((f ree f)) v (3.11) 

It accounts for the direction and length £ of the walls; however, since the symmetry 
£ — > — £ docs not affect the pattern, the orientation of I plays no role. 

The diagonal components M rr and Mgg of this tensor are both of order (£ 2 )- Con- 
versely, the off-diagonal component M r g — Mg r (the tensor is symmetric), and the dif- 
ference M rr — Mgg, together characterise the anisotropy: they are both much smaller 
than (£ 2 )- They vanish when the foam is isotropic; in which case M = ((£ 2 )/2)I is pro- 
portional to /, and only encodes the r.m.s. wall length. Hence the tensor M has two 
strictly positive eigenvalues, also of order (£ 2 ), the larger being in the direction in which 
bubbles elongate. It thus quantifies the average size, as well as the anisotropy direction 
and amplitude, of the network consisting of the bubble vertices (Fig. 2a, top right). In a 
stationary regime, it is constant, as are all averaged physical quantities. 

3.3.2. Reference state 

At any time, the current state of the foam corresponds to a current reference state, 
namely the state that the foam reaches after relaxation (Porte, Berret & Harden 1997). 
This reference state defines the reference texture tensor M . It varies with time; but in 
a stationary regime, it is constant, as are all averaged physical quantities including M. 
The reference state is not necessarily equal to the initial state (Porte, Berret & Harden 
1997). 

In the present experiment, letting the inner wheel rotate back freely is difficult, so 
that in practice we do not have direct access to a reference state. But we do not need 
microscopical details to estimate M : it suffices to make a physical approximation on 
the average properties of the foam at rest. We chose below to assume that the reference 
state is statistically isotropic. We thus approximate Mo by MqI, where the constant 
M = Tr(Ifo) /2 is half the average over all boxes of the squared wall length. We 
implicitly assume that the volume is conserved, which should be correct in Couette flow. 
More precisely, to avoid overestimating Mo due to bubble stretching, we take the average 
between (£} 2 /2 and (£ 2 )/2; since here (£) = 1.40 mm and \J (£ 2 ) = 1.45 mm, we take 
Mo = 1.0 mm 2 . 

3.3.3. Elastic strain tensor 

We can operationally measure the statistical elastic strain tensor U in each box from 
(Aubouy et al. 2003): 

V{f) = l0S ^ (f) - l0S ^. (3.12) 
It obeys (Aubouy et al. 2003) the mathematical requirements for a strain tensor: 
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symmetry properties under translation, rotation and index permutation. It is always 
defined: logM is the tensor with the same eigenvectors as M, but with the logarithm 
of its (strictly positive) eigenvalues. It is the only dimensionless function of M and M 
which is always defined (there is no equivalent of the division for tensors, except when 
they commute). 

It also obeys the physical requirements for an elastic strain tensor. It statistically 
quantifies the deformation rcversibly stored in the present state of the material. Its 
differential dU (and hence the elastic moduli for infinitesimal deformation around the 
current state) depends explicitly only on the current state. 

Moreover, let us consider the case where the (small or large) deformation is elastic and 
affinc. Then: 

e = Te , 

£® £ = £? = T £ 4 7* =f (Y ® £o) T?, 
W = (?® t)=T (l <8> £ ) F* = TTl T*. (3.13) 

If we now assume that M commutes with F, which is the case if for instance M is 

isotropic, then M = M F F . We then prove that U coincides with the true strain 
(Hoger 1987), (Zimmerman 1999): 

, = , = log (M FF ) -log M , 
logM-logMo *\ J 1 #°) 



C^=— = 2 = -\og\FFj=E , (3.14) 

see eq. (3.10). The Finger tensor appears because we base our elastic description on a 
tensor, t ® I = I £*. This constrast with classical elasticity, based on a scalar product 

dx ■ dx = dx 1 dx = dx% F F dxo , where the Cauchy-Green tensor appears (Phan-Thien 
2002). 

A specific feature of the quasistatic regime is that the deformation (and hence the 
elastic strain) remains small, barely visible by eye (Fig. 1), even for large applied shear 
strain (and hence large plastic strain). We thus assume, and check a posteriori (Fig. 7), 

that the components of U are much smaller than one, so we approximate eq. (3.12) by 
its linearisation: 



4. Results 

4.1. Transient regime 

All data in this section are taken from Run 1 (Fig. 1), which has two distinct regimes: a 
transient and a stationary flow (Fig. 3). 

4.1.1. Texture tensor 

In the transient regime, the average wall length does not change much, but the walls 
tend to align with the direction of the rotation (Fig. 1). Hence, it is mostly the cross- 
component M r g which correlates to the shear (Fig. 3). Since 7 = marks the change 
from clockwise to counter-clockwise shear, M r g decreases linearly and changes sign. The 
correlation with shear of both other components, if any, is of order of their fluctuations. 



E. Janiaud and F. Graner 



0.8 



A °- 6 

N 
V 



CM _ . 
<^> 0.4 



0.2 


-0.2 



(a) = 

- h'JP-Cj : - 


: 




(b) 










1 
1 
1 
1 
1 








A 1 

rv.i..vX... i, .£ . •>.•: 








1 v • 









0.5 1 1.5 2 2.5 0.5 1 1.5 2 2.5 

Y Y 



13 




-0.15 t 1 

0.05 0.1 0.15 0.2 0.25 0.3 

y 

Figure 4. Run 1 (transient regime): cross component U r e versus the applied shear strain 7, for 
box 1: closed circles, box 3: closed squares, box 5: open triangles, and box 7: crosses. Since 
7 = marks the reversal of the rotation (and not an initially unstrained state), we expect U r g 
to be zero when 7 = 0.15, within experimental error (± 0.03). A linear fit of U r e versus 7 yields 
the derivative dUre/d'y, plotted in Insert, where the line is a linear fit with zero intercept to 
dUre/d'y versus l/r . 



4.1.2. Elastic strain tensor 

During the transient, no bubble rearrangement (neighbour swapping, also called "Tl 
processes") is visible for 7 < 0.3. This probably corresponds to the elastic regime (Hohlcr, 
Cohen- Addad & Hoballah 1997). The cross-component U r g varies linearly with 7, and 
a linear fit yields its derivative dUre/d'y (Fig. 4). It varies smoothly with r, and a fit to 
r- 2 yields the prefactor (-1600 ± 90) (£) 2 . 

It is interesting to compare dU r g / dj with the expression which would be predicted by 
classical elasticity in a Couette geometry: 

dUre = r 2 r j_ 
c?7 n + ro r 2 

The prefactor of eq. (4.1) is (r 2 r )/[(7i + r )} = (-1680 ±80) {£) 2 , where we have taken 
1.2 mm (the depth of the teeth) as the uncertainty of the wheel radius. This agrees with 
the experimental prefactor of r~ 2 (Fig. 4). 

4.1.3. Yield 

The flow becomes stationary at a value jy which we measure by determining the 
intersection between the transient build-up of M r g fitted by a straight line and the 
stationary plateau fitted by a horizontal line (Fig. 3). It would be interesting to determine 
whether "fy actually marks the onset of irreversibility in the stress-strain relation at 
macroscopic scale (plasticity). 

We can at least observe on the movie (data not shown) the "Tl" bubble rearrange- 
ments, which concern local topology. In the inner boxes, jy correlates with the appear- 
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Figure 5. Run 1 (transient regime): variation of -yy versus box number. Data are from Fig 

(3). The line is a linear fit. 



ance of a large number of Tls. The first isolated Tls appear at values of 7 between 0.26 
and 0.32, of the same order of magnitude as, but significantly smaller than yy. The outer 
boxes lack Tls, except once near 7 ~ 3, well beyond jy and in the stationary regime. 
Though jy is well defined, it marks a saturation in the local shear strain applied to the 
bubbles, which probably oscillates and remains below the actual yield strain. 

The value of "fy increases gradually from 0.35 for the inner box, to 0.55 for the outer 
box (Fig. 3). Within experimental error (± 0.03, corresponding to ± 2 images), this 
variation is linear (Fig. 5). Explaining this variation of jy with r would certainly help 
us to understand the yielding of foams (see section 5.2.2). 

4.2. Stationary regime - average values 

All data of this section are taken from Run 2 (Fig. 2), which allows us to improve our 
statistics, using temporal averages over the stationary flow. 

4.2.1. Texture tensor 

Fig. (6) plots the average M stat of the instantaneous texture tensor M in the stationary 
regime. We can now detect, not only the (small) anisotropy of the foam, M r g, but also 
its small spatial variations. 

M r g decreases across the gap: it is larger in the inner boxes than in the outer ones, 
reflecting the preferred orientation imposed by the shear. M rr decreases, Mgg increases, 
and M rr — Meg changes sign; this sign-change is more difficult to understand (section 
5.1.3). 

4.2.2. Elastic strain 

U r g does not localise near the inner wheel: it varies smoothly with the distance to 
the inner wheel. It decreases as ~ 1/r 2 (Fig. 7), as expected if the foam behaves like a 
continuous elastic medium, see eq. (4.1). 

To compare to the transient regime in Run 1 (which has slightly different bubble and 




Figure 6. Run 2 (stationary flow): texture tensor M stat versus box number. The components of 

M (closed triangles: M rr ; open squares: Mee; closed circles: M r g), again expressed in units of 
(£ 2 ), are averaged over the four quarters of the whole foam, and also over time (2000 pictures). 



wheel diameters), we express the prefactor determined in Fig. (7) in the same units as 
in cq. (4.1): 

(U re )= lel ^^\. (4.2) 
ri + r r z 

The prefactor is 7e ; = 0.171 ± 0.003, see section (5.2.2). 

4.2.3. Deviatoric elastic stress 

In the quasistatic regime, mechanical equilibrium implies that the elastic stress is 
divergence free. Hence we expect that its off-diagonal component, which reduces to the 
capillary contribution, varies as for a continuous elastic medium in a Couette geometry 
(Dcnnin 2004): a c r g P ~ cst/r 2 . Such r -2 dependence occurs experimentally (Fig. 8). 
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Figure 7. Run 2 (stationary flow): tensor U versus box number. The components of U 
(closed triangles: U rr ; open squares: Uee; closed circles: U r e) are derived from the data in 
Fig. (6) through eq. (3.15), using the value Mo = 1.0 mm 2 (see section 3.3.2). The solid line is 
a best fit of U r e to a function est / r 2 . 
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Figure 8. Run 2 (stationary flow): elastic stress tensor a cap versus box number. The deviatoric 
components of !j cap (closed circles: shear stress cr^ p ; open squares: normal stress difference 
arr P — o-gg P ) are derived from eq. (3.5) and expressed in units of A/(£). The solid line is a best 
fit of a c r a v to a function cst/r 2 . 





Figure 9. Run 2 (stationary flow): Hooke relation. The deviatoric part of the capillary 

stress a cap , expressed in units of \/(£), plotted versus the corresponding components of U 
(closed circles: a c °^ versus U r g; open squares: OrV ~ °"ee P versus U rr — Uee). Each data point 
comes from Figs. (7) and (8), i.e. it derives from averages at one position in the foam. The solid 
line is a linear fit through all points, slope 2fj,{£)/\ = 1.17 ± 0.04. 



Again, the spatial variation of a^P — aaa , which changes sign, is more difficult to 
understand (section 5.1.3). 



4.2.4. Non-dimensional shear modulus 

In elasticity, the stress and the deformation are conceptually independent quantities 
(see section 5.1.2). In a continuous medium, they correlate physically; their relation is 
the elastic part of the material's constitutive relation (Landau & Lifschitz 1986). Let us 
focus on the components we measure here, namely the deviatoric terms rr — 99 and r9. 
If the material is linear, the components of strain and stress are proportional to each 
other, the slope being twice the shear modulus (i. If the material is isotropic, rr — 99 and 
r9 play a similar role. In polar coordinates, the Hooke equations becomes: 



aZ P = ^U r6: (4.3) 
(*%>-*%>) = MUtt-Uo6). (4.4) 

Fig. (9) indicates that, in the quasistatic regime, the foam indeed behaves like a contin- 
uous medium: the deviatoric part of the capillary stress cr ca P and of U correlate strongly, 
and details of the microstructure appear only through mesoscopic averages. Moreover, 
the foam is both isotropic and linear, with a shear modulus /i = (0.59 ± 0.02) A/ (£). With 
a smaller or a larger box size, we obtain the same value for /i with a larger error. 
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4.3. Permanent regime - temporal fluctuations 

4.3.1. Time series of M in Run 1 

The stationary flow results from a balance between the increase due to the applied 
shear strain and relaxation due to bubble rearrangements (Langer & Liu 1997). The 
texture tensor and the elastic strain fluctuate around their mean value. 

In Run 1, the short-term fluctuations are much larger near the inner wheel (Fig. 
3a), where shape deformations and bubble rearrangements are both larger (Debregeas, 
Tabuteau & Di Mcglio 2001), than near the outer wheel (Fig. 3b). Fig. (3) docs not 
show any significant physical feature of the fluctuation statistics, e.g. abrupt, large-scale 
relaxations of deformations via correlated, large stress drops. 

4.3.2. Histograms of increments in Run 2 

We perform a more quantitative analysis on Run 2, at a small scale, namely between 
two successive images, corresponding to an increase of d"f in the applied shear strain. 
The increment of the texture tensor is defined as M(j + d-f) — M(-f). The signature 
of "avalanche-like" events at such a small scale would be extremely large fluctuations 
and/or a large asymmetry between positive and negative increments. 

We find no such signature. The increments are almost Gaussian (Fig. lOa-c). The width 
of the histograms decreases by a factor of 10 from inner to outer boxes (Fig. lOd): as we 
mentioned, most large fluctuations occur near the inner wheel. Only M r g displays a slight 
asymmetry between positive and negative increments (Fig. 10b), probably because this 
cross-component is increased by the applied shear strain and decreased by the anisotropy 
relaxation due to Tl events. 

4.3.3. Temporal autocorrelations in Run 2 

The rotation of the inner wheel does not perturb the foam near the outer wheel. In 
the outer boxes, the texture tensor thus decorrelates very slowly. Let us quantify this 
observation using the data of Run 2. 

To analyse the fluctuations at all scales of the shear strain £7, we calculate the temporal 
autocorrelation Cmm of the fluctuations: 

r (M ((M( 7 ) - M gtat ) (M( 7 + Sj) - M stat )) 
Cmm(oj) = 7 -r , (4.5) 



(M( 7 ) - M s1 

where M(j) (resp: M sta t) is the instantaneous value (resp: the average over the stationary 
regime) of the texture tensor M . 

Fig. (11) displays a rapid initial decay in the correlation, followed by a slower decay, 
reflecting the destruction of the initial conditions (Durian 1997), (Lauridscn, Twardos 
& Dcnnin 2002), (Ybcrt & di Meglio 2002), (Pratt & Dennin 2003). An exponential fit 
at small £7 yields the shear strain 7dccorr characterising the initial loss of correlation, 
especially quick for the cross component M r g (Fig. 12): 

M rB = cxp ( . (4.6) 



decorr 



As we expect, the correlation decays much more quickly near the inner wheel, 7d c 
1/30, than near the outer one, 7docorr ~ 1 (insert of Fig. 12). The decorrelation rate 
7decorr decreases exponentially with the distance from the inner wheel. 

The fluctuations of the texture tensor are thus localised near the inner wheel, over 
a distance comparable to 4 bubble diameters. This localisation resembles that of the 
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Figure 10. Run 2 (stationary flow): statistics for the increments of the texture tensor M. 
(a) Histogram of increments of M rr in box 1, again expressed in units of {£ 2 }- Solid line is a 
Gaussian with the same standard deviation and zero average, (b, c) The same for M r g and 

Mgg. (d) Standard deviation of increments of M versus box number (closed triangles: M rr \ 
open squares: Mgg; closed circles: M r g). Since the fluctuations in box 8 are extremely small 
(see Fig. 3), the data are noisy: image filtering and thresholding would have been required to 
obtain a statistically significant value. 



velocity field over 1 bubble diameter (Debregeas, Tabuteau & Di Meglio 2001), and is 
probably the signature of the effect of Tls on the deformation (Dennin 2004). 



5. Discussion 

5.1. Discussion of the strain and stress 

5.1.1. Measurement method of the stress 

Let us discuss the contributions to the stress we do not measure. 
The dissipation is probably dominated by the friction on both glass plates, which we 
expect to be much larger than the internal viscous stress (Cantat, Kern & Delannay 
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Figure 11. Run 2 (stationary flow): temporal autocorrelation Cmm of texture tensor M (eq. 
4.5), versus S'y, in boxes 1, 3, 5, 7 (closed triangles: M rr ; open squares: Mee; closed circles: 
M r g). We recall that the total applied shear 7 is proportional to time, not to the local shear in 
each box (eq. 2.1). 



2004). Whatever its physical origin, it should on average balance the input power. It is 
thus never negligible, especially during the relaxation after each Tl, see section (2.1). 
However, the results we present here, which concern only the elastic terms, do not require 
us to measure the dissipation. 

We measure the deviatoric elastic stress using eq. (3.5), which is formally equivalent 
to the stress in a network where sites interact through a two-body force. Foams do 
not include only two-vertex interactions. Their stability is due to the pressure inside 
each bubble, which relates to the bubble volume. Hence pressure is a function of all 
wall curvatures and all the vertex positions. It is an n-body interaction, which we can 
decompose into two-point interactions in principle (Alexander 1998), but not really in 
practice. 

It is not important here to measure the pressure inside bubbles, for the following 
reasons (Graner 2002). Since the pressure stress is isotropic, it does not contribute to the 
deviatoric elastic stress. Hence the shear modulus is independent of the pressure field. 
Finally, the compression modulus of a foam is typically that of an ideal gas, which scales 




5y 

Figure 12. Run 2 (stationary flow): autocorrelation C l M r9 M I . e of the cross-component M r e. Box 
1: crossed boxes, box 2: open circles, box 3: closed squares, box 4: open triangles, box 5: 
closed triangles, box 6: closed circles, box 7: open squares, box 8: diamonds. Insert: initial 
decorrelation rate 7,7 C c 0rr for each box (eq. 4.6). 



like the pressure inside bubbles ~ 10 5 Pa (multiplied by h, when expressed as a 2D 
compression modulus) . This modulus is not very interesting and is much larger than the 
shear modulus of the foam, which is thus effectively incompressible (Weaire & Hutzler 
1999). 

The line tension is not known, but unimportant for our present purposes. The stress, 
and hence the shear modulus, are expressed in non-dimensional form, as required for 
comparison to other foams. A direct measurement of A can complement the local image 
measurement when necessary, for instance to predict the actual force exerted by the 
flowing foam on an obstacle, or to validate eq. (3.5) by checking its agreement with an 
independent, macroscopic measurement of the force (Courty et al. 2003). 

5.1.2. Difference between a cap and M 

Both W cap and M are measured from the same image, and their definitions (eqs. 3.5, 
3.11) seem similar. However, they are really independent both physically and mathemat- 
ically. 

Physically, their significance and status are different. The value of cf cap can be deter- 
mined from image analysis only if the sites interact through a two-body force r, and if 
we explicitly know the form of this force f{i) with respect to the distance between sites. 
The value of a cap is sensitive to physical forces: first, it doubles if the density of links 
doubles; second, it does not change if we cut a link in two by artificially introducing a 
new site at its middle; third, in foams, it formally depends on wall curvatures. 

On the other hand, M is purely geometric. It is a characteristic of the pattern, in- 
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dependent of physical interactions; it can be denned and measured for various kinds of 
cellular patterns for which no f can be defined, including solid foams, and even grain 
boundaries in crystals (Durand, Weiss & Graner 2004). The value of M is sensitive to 
the detailed network topology: first, it does not change if the density of walls doubles; 
second, it decreases if we cut a wall in two by artificially introducing a new vertex at 
its middle; third, in foams, it is insensitive to the (in-plane or out-of-plane) curvature of 
walls. 

Mathematically, 7j cap and M are independent in all cases, except when f oc £, i.e. 
a harmonic spring, with a length I much greater than its length at equilibrium. For 
instance, writing the stress of an isolated polymer molecule in a shear solvant as t ® I 
is legitimate (Bird, Armstrong & Hassager 1987). However, foams, for which t = const, 

fall in the general case where 7j cap and M are mathematically different averages. The 
difference is even more visible in a 3D foam, where cf cap is an integral over a surface 
(Batchelor 1970), while M remains based on vectors (Aubouy et al. 2003). 

Moreover, M does not even need be defined on the vectors I linking neighbour vertices. 
Definitions other than eq. (3.11) could be acceptable in principle (Dollct et al. 2004). For 
instance, we could use the vector g which links the centers of mass of neighbouring 
bubbles, and construct Mas (/ g<%>g). This would apply to bubble rafts (Dennin 2004) , 
granular materials (Kruyt 2003), or 3D foams (Kraynik, Reinelt & van Swol 2003). This 
would change both M and Mo by a multiplicative constant: the elastic strain U (eqs. 
3.12, 3.15), and hence the shear modulus fi (eq. 4.4) would almost not change. 

We emphasise that a cap and U are also independent. They represent different aver- 
ages of microscopic details, so we cannot express one as a function of the other (Kruyt 
2003), (Goldhirsch & Goldcnbcrg 2002), except empirically, through numerical simula- 
tions or physical approximations approximations (Kruyt 2003). The difference is even 
more visible at large deformations: the logarithmic term of eq. (3.12) deviates from the 
linear approximation when the maximum value of the elastic shear strain (0.15 in Fig. 7) 
reaches large values, e.g. 0.4 (Courty et al. 2003) or 0.6 (Asipauskas et al. 2003). There is 
no trivial relation between them, except in a regular pattern, where averages are replaced 
by exact identities; for instance, in the case of a honeycomb, their relation is exactly as 
expected (Dollet et al. 2004). 

5.1.3. Sign of normal differences 

The normal differences M rr — Mgg and cr™ p — cr^ p change sign in the middle of the 
gap, between boxes 4 and 5 (Figs. 6, 7, 8). 

Normal stresses appear in elastic solid (Poynting 1909) or granular (Bagnold 1941) 
materials sheared at constant volume (this is equivalent to "dilatancy" : volume increase 
under shear at constant pressure (Weaire & Hutzler 2003)). In complex fluids such as 
foams, these normal stresses exist but are difficult to measure (Khan, Schnepper & Arm- 
strong 1988) and to predict (Weaire & Hutzler 2003). Predictions and 2D simulations of 
their average are accurate: under steady shear, the pressure of a foam with fluid fraction 
$ = 5% increases by 15% (Weaire & Hutzler 2003). This is the order of magnitude of the 
spatial variations of M rr and Mgg (Fig. 6). We can not predict their anisotropy (hence 
the normal stress differences), but normal stresses due to dilatancy are a second-order 
effect in the strain (Larson 1997): they are likely smaller that the large value of M rr — Mgg 
we observe, and might be masked by dynamical effects due to bubble rearrangements. 
But we can imagine other explanations. 

First, this change could reflect an initial anisotropy resulting from the procedure used 
to fill the cell with foam (Debregeas, Tabuteau & Di Meglio 2001), which the preparatory 
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rotation does not suffice to relax in the outer boxes (Debregeas, personal communica- 
tion). The initial state is not necessarily a zero-stress state. More generally, trapped 
stresses make notoriously difficult to measure normal stress differences (Khan, Schnep- 
per & Armstrong 1988), and carefully relaxed simulations are more reliable (Kraynik, 
Reinelt & van Swol 2003). The present measurement is probably one of the first direct 
experimental ones, along with that of Ref. (Labiausse, Hohler & Cohen- Addad 2004). 

Second, it could be an artefact due to our measurement method, which docs not 
accounts for the out-of-plane curvature of bubble walls (in the third dimension). Since 
the foam is not perfectly dry, the sides of a wall can have slightly different curvatures, 
so that their contributions would not exactly cancel out. 

Third, this change could arise from a deviation from axisymmetric Couette flow, due 
to a secondary recirculation with radial velocity components. A signature of the discrete 
and disordered nature of the material is the appearance of solid body rotation of bubble 
clusters. It results in vortices, generic for amorphous elastic materials, generated by the 
"backflow" of the non-affine displacements (Tanguy et al. 2002). Their lifetime is as short 
as the time the internal wheel takes to move one bubble diameter (Debregeas, Tabuteau 
& Di Meglio 2001). 

5.2. Discussion of the stationary flow 
The present measurements clarify a few open debates. 

5.2.1. To localise or not to localise? 

The experimental set-up, the fluid fraction, the friction on the plates, and the boundary 
conditions at the wheel, are those of Ref. (Debregeas, Tabuteau & Di Meglio 2001). How 
can the same data alternatively display and not display localisation? 

There is no fundamental incompatibility: different fields have different spatial varia- 
tions. Ref. (Debregeas, Tabuteau & Di Meglio 2001) measures the velocity field, and the 
Tl rate; while we measure the elastic stress and strain fields. In a quasistatic regime, 
where the foam is almost always close to mechanical equilibrium, we expect the stress to 
vary smoothly across the gap, as observed here and in simulations (Kabla & Debregeas 
2003). The elastic strain, which correlates to the elastic stress, does not localise, although 
the plastic strain, which correlates to the velocity gradient, does. 

5.2.2. Yield 

The value of the yield strain and stress, characterised e.g. by the equality of the real 
and imaginary parts of the complex shear modulus (Hohler, Cohcn-Addad & Asnacios 
1999), is not completely understood yet. It generally depends on the shear rate (Coussot 
et al. 2002), (Rouyer et al. 2003), but not in the present quasistatic regime. 

We expect that, to store a deformation 7 e ;, requires us to apply an equivalent shear 
strain yy. Since Run 1 corresponds to an inversion from counter-clockwise to clockwise 
rotation, we expect that yy = ^lei- We observe an excellent agreement in box 1, in which 
yy = 0.35 ± 0.02 (Fig. 3) and 7eJ = 0.171 ± 0.003 (Fig. (7) and cq. (4.2)). But yy ^ 2 7e , 
in other boxes, for larger yy (Fig. 5). Thus the deformation stored in the foam changes 
with 7 even above yy. Kabla and Debregeas reach the same conclusion: their simulations 
indicate that, after the foam yields, Tls localise in the shear band only after a certain 
delay time (Kabla & Debregeas 2003). 

5.2.3. Averages and fluctuations 

Our data show that stress and strain have a well defined, physically meaningful average 
at the scale of a representative volume element, while their fluctuations remain small. A 
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coarse-grained shear modulus can be denned and measured, unlike e.g. in granular media, 
where the importance of fluctuations at large scales is still an open debate (Goldhirsch 
& Goldenberg 2004). 

More precisely, the statistics of the texture tensor are nearly Gaussian. Visual inspec- 
tion of the images confirms that Tls remain isolated. They correlate (Kabla & Debregeas 
2003), but over a very small space and time range, probably due to the experimental 
conditions: disordered foam, quasistatic regime, large dissipation due to top and bot- 
tom plates. We do not observe extremely large fluctuations, nor sudden relaxations of 
deformations and stress. 

Averages seem to dominate the physics Fig. (3), but fluctuations too are important. 
For instance, the width of the histogram of the velocity or stress field fluctuations might 
constrain the velocity field, and explain why shear-banding occurs preferentially near 
the inner wheel, due to a feedback loop on which fluctuations grow (Kabla & Debregeas 
2003). 

Both the average and fluctuations of Tl rate manifest themselves in the texture tensor. 
The signature of the average Tl rate is the existence of a stationary regime, where the 
decrease in deformation due to Tls balances the increase due to applied shear. The 
signature of Tl rate fluctuations is the exponential localisation of the texture tensor 
fluctuations. 



5.2.4. Measurement of the shear modulus 

Fig. (9) suggests a few remarks. First, the classical method to measure a shear mod- 
ulus consists in increasingly deforming a sample, and measure the stress and strain in 
successive states. We use here another method (Asipauskas et al. 2003). Each point 
in Fig. (9) comes from a different region of the same foam. We thus observe different 
deformation states simultaneously from a single heterogeneous image. Averaging over 
successive images reduces noise (Asipauskas et al. 2003), but is not necessary (Courty et 
al. 2003), (Durand, Weiss & Graner 2004). 

Second, when determining the constitutive relation for plastic materials, the elastic 
stress is classically plotted versus the total applied shear strain: see Refs. (Jiang et al. 
1999), (Reinelt & Kraynik 2000), (Weaire & Hutzler 2003) to quote a few. Different ap- 
plied strains correspond to the same stress. Here we observe a one to one correspondance 
between elastic stress and elastic strain (Fig. 9). The foam is elastic, although the plastic 
strain increases steadily, and the foam is well beyond the yield strain. 

Third, although we do not know the value of A, we can compare our non-dimensional 
measurement of shear modulus /x with published values. An ideal two-dimensional dry 
monodisperse foam, where bubbles form a regular honeycomb lattice (Princen 1983), 
(Khan & Armstrong 1986), (Weaire & Hutzler 1999), has a shear modulus fj,h c = ^ = 

0.577-^-. Our \i is slightly larger, probably due to disorder; here )jp = 1.08. The non- 
zero fluid fraction tends to decrease /i (Princen 1983), (Khan & Armstrong 1986), (Ma- 
son, Bibette, & Weitz 1995), (Weaire & Hutzler 2003). In a non-quasistatic, and more 
heterogeneous, stationary flow, either through a constriction or around an obstacle, a 
dry foam (fluid fraction much lower than 5%), even with a small bubble size dispersity 
< 5%, displays a shear modulus 20% higher than /z^ c , due to the stretched bubble shape 
(|0 = 1.22 (Asipauskas et al. 2003). 
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6. Conclusion 

To summarise, we re-analyse the data generously provided by Georges Debregeas on a 
foam sheared in a two-dimensional Couette geometry, with the velocity field localised near 
the inner wheel (Debregeas, Tabuteau & Di Meglio 2001). Since the flow is quasistatic, 
the deformation of the bubbles is very small. We take advantage of the foam's invariance 
under rotation (and, in the stationary regime, of its time invariance) to average over many 
bubbles. We then measure the texture tensor of the foam, a general tool to quantify the 
deformation of 2D or 3D networks and cellular patterns (Aubouy et al. 2003). 

For small applied shear strain, a transient regime has the cross-component of the 
texture tensor increase linearly with shear, with no bubble rearrangements ("Tl"), and 
the foam seems to behave like a classical elastic, continuous medium. The only large 
scale manifestation of the discrete, cellular nature of the foam is that, in a mixture of 
two fluids (soap solution and air), the surface tension of bubble walls produces elastic 
behaviour and a shear modulus. 

The deformation saturates first near the inner wheel, then gradually in more distant 
boxes, at an applied shear strain jy- This onset of plastic behaviour correlates, but does 
not strictly coincide with the appearance of Tls. 

The foam then enters a stationary regime, where the deformation increase due to shear 
and the relaxation due to Tls balance each other. The fluctuations of the texture tensor 
around its mean value remain small. Despite the Tls, the elastic behaviour of the foam 
again resembles that of a classical elastic, continuous medium. The maximum shear strain 
that the foam can sustain in the quasistatic regime is "f e i, compatible with the value of 
7y. The strain and stress are highly correlated, even proportional to each other; the 
shear modulus is slightly higher than for a honeycomb with similar average bubble wall 
length. The deviatoric components of the elastic strain and stress tensors, as well as the 
texture tensor, do not localise. The only trace of localisation is the much quicker loss of 
correlation of the texture tensor fluctuations near the inner wheel. 

Future studies aim at understanding the origin of the value of the yield strain, the 
change of sign in the normal stress and strain differences, the precise role of boundary 
conditions, and the relation to granular materials. They include the measures of pressure 
and dissipation, and the analytical calculation of the shear modulus. 

Acknowledgments: We are extremely grateful to G. Debregeas for providing us with 
many recordings of his published and unpublished experiments, and for detailed discus- 
sions. We would also like to thank M. Dennin, R. Holder, A. Kabla, A. Kraynik, N. 
Kruyt, J. Scheibert for discussions, M. Aubouy, B. Dollct, F. Elias and J. A. Glazier for 
critical reading of the manuscript, O. Cardoso for help with image analysis. This work 
was partially supported by CNRS ATIP 0693. 

REFERENCES 
Alexander S. 1998, Phys. Reports, 296, 65. 

Asipauskas M., Aubouy M., Jiang Y., Graner. F. & Glazier J. A. 2003, Granular Matt, 
5, 71. 

Aubouy M., Jiang Y., Glazier J. A. & Graner F. 2003, Granular Matt., 5, 67. 
Bagnold R.A. 1941, The Physics of Blown Sand and Desert Dunes, Methuen, London. 
Ball R. & Blumenfeld R. 2002, Phys. Rev. Lett, 88, 115505. 
Batchelor G.K. 1970, J. Fluid Mech., 41, 545. 

Bird R.B., Armstrong R.C. & Hassager O. 1987, Dynamics of polymeric liquids, John 

Wiley & Sons, NY (USA). 
Blumenfeld R. 2004, Physica A, 336, 361. 

Cantat I., Kern N. & Delannay R. 2004, Europhys. Lett. 65, 726. 



26 



E. Janiaud and F. Graner 



Courty S., Dollet B., Elias F., Heinig P. & Graner F. 2003, Europhys. Lett. 64, 709. 
Coussot P., Raynaud J.-S., Bertrand F., Moucheront P., Guilbaud J. -P., Huynh H.T., 

Jarny S. & Lesueur D. 2002, Phys. Rev. Lett. 88, 218301. 
Cox S.J., Weaire D., & Vaz M.F. 2002, Eur. Phys. J. E, 7, 311. 
Debregeas G., Tabuteau H. & Di Meglio J.-M. 2001, Phys. Rev. Lett, 87, 178305. 
Dennin M. 2004, 

Dollet B., Mecke K., Aubouy A. and Graner F. 2004, in preparation. 
Durand G., Weiss J. & Graner F. 2004, Europhys. Lett, 67, 1038. 
Durian D. J. 1997, Phys. Rev. E, 55, 1739. 

Elias F., Flament C, Glazier J. A., Graner F. & Jiang Y. 1999, Philos. Mag. B, 79, 
729. 

Farahani K. & Naghdabadi R. 2000, Int. J. Solids Struct., 37, 5247. 

Geminard J.-C, Zywocinski A., Caillier F. & Oswald P. 2004, Phil. Mag. Lett. 84, 199. 
Goldhirsch I. & Goldenberg C. 2002, Eur. Phys. J. E, 9, 245. 
Goldhirsch I. & Goldenberg C. 2004, Granular Matter, 6, 87. 
Graner F., Jiang Y., Janiaud E. & Flament C. 2001, Phys. Rev. E, 63, 011402. 
Graner F. 2002, Two-dimensional fluid foams at equilibrium, in Ref. (Mecke & Stoyan 2002), 
pp. 187-214. 

Hoger A. 1987, Int. J. Solids Struct, 23, 1645. 

Hohler R., Cohen-Addad S. & Asnacios A. 1999, Europhys. Lett., 48 , 93. 

Hohler R., Cohen-Addad S. & Hoballah H. 1997, Phys. Rev. Lett., 79, 1154. 

Jiang Y., Asipauskas M., Glazier J. A., Aubouy M. & Graner F. 2000. In Eurofoam 2000 

ed. Zitha P., Banhart J. & Verbist G., p. 297, MIT Verlag, Bremen. 
Jiang Y., Swart P.J., Saxena A., Asipauskas M. & Glazier J. A. 1999, Phys. Rev. E, 59, 

5819. 

Kabla A. & Debregeas G. 2003, Phys. Rev. Lett., 90, 258303. 

Kern N. & Weaire D. 2003, Phil. Mag. 83, 2973. 

Khan S.A. & Armstrong R.C. 1986, J. Non-Newt. Fluid Mech., 22, 1. 

Khan S.A. & Armstrong R.C. 1987, J. Non-Newt. Fluid Mech., 25, 61. 

Khan S.A., Schnepper C.A. & Armstrong R.C. 1988, J. Rheol, 32, 69. 

Kraynik A., Reinelt D., & van Swol F. 2003, Phys. Rev. E, 67, 031403. 

Kruyt N.P. 2003, Int. J. Solids Struct, 40, 511. 

Kruyt N.P. & Rothenburg L. 1996, J. Appl. Mech., 118, 706. 

Kruyt N.P. & Rothenburg L. 2002, Int. J. Solids Struct, 39, 311. 

Labiausse V., Hohler R. & Cohen-Addad S. 2004 in Proceedings of Eufoam 2004, Adler 

M. ed., Elsevier, to appear. 
Landau L. D. & Lifschitz E. M. 1986, Theory of Elasticity, Reed, Oxford, 3 rd ed.. 
Danger S.A. & Liu A.J. 1997, J. Phys. Chem. B 101, 8667. 
Larson R. 1997, J. Rheol. 41, 365. 

Lauridsen J., Twardos M. & Dennin M. 2002, Phys. Rev. Lett, 89, 098303. 
Liao C.L., Chang T.P., Young D.H. & Chang C.S. 1997, Int. J. Solids Struct, 34, 4087. 
MACOSKO C. 1994, Rheology: principles, measurements and applications, Wiley- VCH, New- 
York. 

Mason T., Bibette J., & Weitz D. 1995 Phys. Rev. Lett. , 75, 2051. 

Mecke K. & Stoyan D. eds. 2002, Morphology of Condensed Matter - Physics and Geom- 
etry of Spatially Complex Systems, Lecture Notes in Physics 600, Springer, (Heidelberg), 
proceedings of the 2nd conference "Spatial Statistics and Statistical Physics" , Wuppertal 
(Germany), March 2001. 

Phan-Thien N. 2002, Understanding viscoelasticity, Springer, Berlin. 

Porte G., Berret J.F. & Harden J.L. 1997, J. Phys. II France, 7, 459. 

Poynting J.H. 1909, Proc. Roy. Soc. London, A82, 546. 

Pratt E. & Dennin M. 2003, Phys. Rev. E, 67, 051402. 

Princen H. M. 1983, J. Coll. Interf. Sci. 91 160. 

Reinelt D. A. & Kraynik A. 2000, J. Rheology, 44, 453. 

Rosenkilde C.E. 1967, J. Math. Phys., 8, 84. 



Local measurement of bubble deformation 



27 



Rouyer F., Cohen-Addad S., Vignes-Adler M. & Hohler R. 2003, Phys. Rev. E, 67, 
021405. 

Scheibert J., Debregeas G. 2004, in preparation. 

Tanguy A., Wittmer J. -P., Leonforte F. & Barrat J.-L. 2002, Phys. Rev. B, 66, 174205. 
Tanner R.I. & Tanner E. 2003, Rheo. Acta, 42, 93. 

Weaire D. & Hutzler S. 1999, Physics of Foams, Oxford Univ. Press, Oxford. 
Weaire D. & Hutzler S. 2003, Phil. Mag., 83, 2747. 

Weaire D., Kern N., Cox S.J., Sullivan J.M. & Morgan F. 2004, Proc. Roy. Soc: Math. 

Phys., 460, 569. 
Ybert C. & di Meglio J.-M. 2002, C.R. Physique, 3, 555. 

Zimmerman J. A. 1999, PhD thesis, Univ. Stanford, chap. 4, pp. 67-106, unpublished. 



